Survival Probabilities in Coherent Exciton Transfer with Trapping 
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In the quest for signatures of coherent transport we consider exciton trapping in the continuous-time quantum 
walk framework. The survival probability displays different decay domains, related to distinct regions of the 
spectrum of the Hamiltonian. For linear systems and at intermediate times the decay obeys a power-law, in 
contrast to the corresponding exponential decay found in incoherent continuous-time random walk situations. 
To differentiate between the coherent and incoherent mechanisms, we present an experimental protocol based 
on a frozen Rydberg gas structured by optical dipole traps. 
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Recent years have seen an upsurge of interest in coherent 
energy transfer, given the experimental advances in manip- 
ulating and controlling quantum mechanical systems. From 
the theoretical side, such investigations are of long standing; 
see, e.g., [1]. Here, tight-binding models, which model coher- 
ent exciton transfer, are closely related to the quantum walks 
(QW). As their classical random walk (RW}^ counterpart, QW 
appear in two variants: discrete-time QW [2] and continuous- 
time QW (CTQW) 13J. Experimental implementations have 
only recently been proposed for both QW variants, based, e.g., 
on microwave cavities I'T], ground state atoms |5] or Rydberg 
atoms in optical lattices, or the orbital angular momentum 
of photons 101 ■ 

An appropriate means to monitor transport is to follow the 
decay of the excitation due to trapping. The long time de- 
cay of chains with traps is a well studied problem for clas- 
sical systems @|: for an ensemble of chains of different 
length with traps at both ends the averaged exciton survival 
probability has a stretched exponential form exp(— M'*'), with 
A = 1/3 (see, e.g., 19]). In contrast, quantum mechanical 
tight-binding models lead to A = 1/4 llOllllll . However, up 
to now only little is known about the decay of the quantum 
mechanical survival probability at experimentally relevant in- 
termediate times. 

Here we evaluate and compare the intermediate-time de- 
cays due to trapping for both RW and QW situations by 
employing the similarity of the CTRW and the CTQW for- 
malisms. Without traps, the coherent dynamics of excitons on 
a graph of connected nodes is modeled by the CTQW, which 
is obtained by identifying the Hamiltonian Hq of the system 
with the CTRW transfer matrix Tq, i.e., Hq = — Tq; see 
e.g. [3, 12] (we will set = 1 in the following). For undi- 
rected graphs. To is related to the connectivity matrix Aq of 
the graph by Tq = — Aq, where (for simplicity) all transmis- 
sion rates are taken to be equal. Thus, in the following we take 
Ho = Ao. The matrix Ao has as non-diagonal elements A^^^ 
the values — 1 if nodes k and j of the graph are connected by 
a bond and otherwise. The diagonal elements of Ao 
equal the number of bonds fj which exit from node j. By 
fixing the coupling strength between two connected nodes to 
I ^ '■^^ scale is given in units of [fi,/i/^°j]. For 



the Rydberg gases considered in the following, the coupling 
strength is roughly /^^ ?J 1 MHz, i.e., the time unit for 
transfer between two nodes is of the order of a few hundred 
nanoseconds. 

The states |j) associated with excitons localized at the 
nodes j {j — 1, . . . , N) form a complete, orthonormal ba- 
sis set (COBS) of the whole accessible Hilbert space, i.e.. 



{k\j) 



and '^i^\k ){ k\ = 1. In general, the time 



evolution of a state |j) starting at time io = is given 
by = exp(— iHo<)|j); hence the transition amplitudes 
and the probabilities read afcj (i) = (fc| exp(— iHot)|j) and 
TTkj{t) = I afcj (i) I ^, respectively. In the corresponding classi- 
cal CTRW case the transition probabilities follow from a mas- 
ter equation as (t) = (fc| cxp(Toi)lj) 12,12]. 



Consider now that out of the N nodes M are traps with 
M < N; we denote them by m, so that m E A4, with A4 C 
{1, . . . , N}. We incorporate trapping into the CTQW for- 
malism phenomenologically by following an approach based 



g an approacn i 

on time dependent perturbation theory lliol [ill Q. The 
new Hamiltonian is H = Ho + iT, where the trapping 
operator iT has at the trap nodes m purely imaginary di- 
agonal elements ir,„,„, which we assume to be equal for 
all m (Fmm = r > 0), and is zero otherwise. As a 
result, H is non-hermitian and has N complex eigenval- 
ues. El = e; — iji (I = 1,...,N). In general, H has 
N left and N right eigenstates 1$;} and ($;|, respectively. 
For most physically interesting cases the eigenstates can 
be taken as biorthonormal, = Sw, and complete, 

J2iLi = 1; see, e.g., Ref. flT]. Moreover, we have 

= {^i\k). Thus, the transition amplitudes can be cal- 
culated as akj{t) = X); exp[-7/t] exp[-ie/t](fc|$;)(<l;|j); 
here the imaginary parts 7; of Ei determine the temporal de- 
cay of TTkjit) = \akj{t)\^- 

In an ideal experiment one would excite exactly one node, 
say j ^ A4, and read out the outcome iTkj (t), i.e., the proba- 
bility to be at node k ^ Ai at time t. However, it is easier to 
keep track of the total outcome at all nodes k ^ A4, namely, 
J2k(^M '^kj{t)- Since the states |fc) form a COBS we have 
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T,k^M = 1 - EmeA4 \m){m\, which leads to: 



N 



N 



-i(Ei-E;,)t 



1=1 
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Ll' = l 



^ {j\^i,){^i,\m){m\<^i){^i\j). (1) 



By averaging over all j ^ M., the mean survival probability 

Ha/ (t) = Ej^M Efc0A4 T'kj {t) is given by 



1 



N 



1 



N - M 



N 



1=1 



meM 



(2) 



For CTRW we include trapping in a formally similar fash- 
ion as for the CTQW. Here, however, the classical transfer 
matrix Tg is modified by the trapping matrix T, such that the 
new transfer matrix is T = To — F, [15]. For a single linear 
system with traps at each end, the mean survival probability 
PM{t) = whiT.jt^MT.ki^MPkAi) decays exponentially 
at intermediate and at long times [15.1. As we proceed to show, 
the decays of njv/(i) and PM{t) are very different, thus allow- 
ing to distinguish experimentally whether the exciton transfer 
is coherent or not. 

For long t and small M/N , Eq. dU simplifies considerably: 
At long t the oscillating term on the right hand side drops out 
and for small M/N we have 2 I]meM(*^'l™)("^l*^') ^ ^■ 
Thus, IImIO is mainly a sum of exponentially decaying 
terms: 



1 



N 



N ~ M 



^exp[-27/i]. 



(3) 
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Asymptotically, Eq. ^ is dominated by the 7/ values clos- 
est to zero. If the smallest one, 7min^ is well separated from 
the other values, one is led for i ^ 1/ 7inin to the exponential 
decay found in earlier works, IImC^) = exp(— 27inini) lflll[ . 

Such long times are not of much experimental relevance 
(see also below), since most measurements highlight shorter 
times, in which many 7; contribute. In the corresponding en- 
ergy range the 7/ often scale, as we show in the following, so 
that in a large I range 7/ ^ aZ'*. The prefactor a depends only 
on r and N 111 l|[ . For densely distributed 7; and at intermedi- 
ate times one has, from Eq. (O, 



^M{t) ~ I dx e 



dy 



(2ai)i/p 



(4) 



The envisaged experimental setup consists of clouds of 
ultra-cold Rydberg atoms assembled in a chain over which an 
exciton migrates; the trapping of the exciton occurs at the ends 
of the chain. The dipolar interactions between Rydberg atoms 
depend on the mutual distance R between the nodes as R^^. 



Now, CTRW over a chain of regularly arranged sites lead both 
for nearest-neighbor steps and for step distributions depending 
on R as R^'^, with 7 > 3, to a standard diffusive behavior and, 
therefore, belong to the same universality class, see e.g. il6ll . 
The reason is that in one dimension for 7 > 3 the first two 
moments, (i?) and (i?^), are finite. Thus, although the quan- 
titative results will differ, the qualitative behavior is similar. 
Hence, we focus on a nearest-neighbor tight-binding model 
and consider a chain of length N with two traps (Af — 2) 
located at its ends (to = 1 and m = N), Q. The CTQW 
Hamiltonian thus reads 

N 

H = J2(2\n){n\-\n-l){n\-\n + l){n\) 



n=l 

+iT |TO)(m| 

m=l,N 



(5) 
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FIG. 1: (Color online) Imaginary parts 7; (dots) in ascending order 
for — 100 and F = 1. Note the shortened y axis. The inset shows 
7i in log-log scale for / = 10, . . . , 90. 

In Fig. [T]we show the spectrum of 7; for = 100 and 
F = 1; the double logarithmic plot (see inset) demonstrates 
that scaling holds for 10 < Z < 60, with an exponent of 
about ^ = 1.865. In this domain 7; G [0.0012,0.012], 
which translates to experimentally accessible coherence times 
of about 10 — 100/is. For comparison the smallest decay rate 
is 7min = 7.94 X 10^^, which corresponds to experimentally 
unrealistic coherence times of the order of tenths of seconds. 

The corresponding transfer matrix of the classical CTRW 
reads 

N 

T = -^(2|n)(n|-|n-l)(n|-|n + l)(n|) 

n=l 

-r J2 (6) 

7n=l,N 

In Fig.|2]we compare the classical PM{t) to the quantum me- 
chanical survival probability 11^/ (t) for a linear system with 
N = 100 and F = 1. Evidently, P/i/(t) and Ha/ (t) differ 
strongly: the PM{t) decay established for CTRW is prac- 
tically exponential. IIa/ (t), on the other hand, shows two 
regimes: a power-law decay at intermediate times (upper red 
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FIG. 2: (Color online) Temporal decay of IlM{t) (solid black lines) 
and PM{i) (short dashed green lines) for A'^ = 100 and F = 1 
in double logarithmic scales (upper three curves) and in logarithmic 
scales (lower three curves). Indicated are the fits to IlM{t) (long 
dashed lines) in the intermediate (upper red) and the long (lower 
blue) time regime. 



curve) and an exponential decay (lower blue curve) at very 
long times. 

We now turn to the parameter dependences of IIm {t). Fig- 
ure |3] displays the dependence of IiM{t) on N . We note that 
the scaling regime, where IiM{t) t^^^^, gets larger with 
increasing N. The cross-over to this scaling region from the 
domain of short times occurs around t w N/2. For larger N 
and in the intermediate time domain IIm {t) scales nicely with 
N. In this case, the power-law approximation [Eq. (|4|i] holds 
and by rescaling / to / /N we get from Eq. Q that 



TlAiit) 



I 



(7) 

where we used that a ^ N^'^ for a linear system llllll . Thus 
when rescaling Itol/N, time has to be rescaled by the factor 
jV-(3-M). Indeed, all curves where a power-law behavior can 
be justified fall on a master curve; see the inset in Fig. [3] 
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FIG. 3: (Color online) A^-dependence of \lM(t) for F = 1; A'^ in- 
creases in steps of 10 from 20 (blue line) to 100 (green line). The 
inset shows IlM{t) versus the rescaled time t/N^^''''. 

The temporal decay does not only depend on N but also on 
r. Figure g] shows HmIO for iV = 50 and different P. For 
values r ^ 1 (green lines) and F <C 1 (black lines) the curves 



shift to longer times. Values of F close to 1 (blue lines) lead 
to the quickest decay. Note that these values are of the same 
order as the coupling strength between the non trap nodes, 

Hjj±i = -1. 
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FIG. 4: (Color online) F-dependence of TlM{t) for intermediate t 
and iV = 50. 

An experimental implementation of the described system 
has to meet several criteria. A single node must represent 
a well-defined two-level system to ensure coherent energy 
transfer while at the same time a mechanism is needed to 
trap an exciton with a controllable trapping efficiency. Fur- 
thermore, the chain must be static with negligible motion and 
should allow for spatially selective excitation and detection of 
the exciton. These demands rule out many possible candidates 
for an experimental realization of CTQW. A frozen Rydberg 
gas [17] can meet all of the above demands by combining 
the rich internal structure of highly excited atoms whith the 
full quantum control over the external degrees of freedom that 
is available in up-to-date experiments with ultracold atoms. 
The internal structure of Rydberg atoms provides both decou- 
pled two-level subsystems and tunable traps, while the pro- 
nounced Stark shift allows to selectively address single sites in 
a chain when an electric field gradient is applied. At the same 
time, experimentally accessible temperatures below 1 /iK en- 
sure that the thermal motion is negligible. 

Our scheme starts from a cloud of laser-cooled ground state 
atoms prepared in a chain of optical dipole traps [18]. Each 
site represents one node with distances between sites of 5 to 
20 /im. For an experimentally achievable extension of 1 mm 
this translates into approximately 100 nodes. All nodes are ex- 
cited to Rydberg states exploiting the dipole blockade mech- 
anism to ensure a single Rydberg excitation per node fl9ll 
which avoids many -body effects [20]. A two-level system is 
realized by a transition between energetically isolated states, 
i.e., by low-angular-momentum states which exhibit a large 
quantum defect, e.g., nS ^ nP . A number of ex- 
periments has revealed the coherent character of this pro- 
cess ll2dh . By contrast to low-^ states, states with angular 
momenta £ > 3 have no quantum defect and are degenerate. 
This allows to construct an exciton trap with the transitions 

n'D ^ n"F ^ n"£{e > 3), where the first transition is the 
dipole transition providing the coupling to neighboring nodes 
II25II while the second transition, driven by a radio-frequency 
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(rf) field, represents the trap and decouples this site from the 
energy transfer, as the large degeneracy of the high-£ states en- 
sures an efficient suppression of the coupling back to the n"F 
state 1 26]. By changing the strength of the driving rf field, 
the trapping efficiency can be tuned. The population of the 
n"£ state is directly proportional to 1 — 11^(1) and can be 
determined by state selective field ionization 12 111 . In an ex- 
periment the central nodes would be prepared in the S state 
and the trap nodes in the D state. A single S node is swapped 
to P through a microwave transition in an electric field gradi- 
ent which makes the resonance S^P position sensitive. This 
is equivalent to exciting a single exciton. The energy transport 
is started by removing the field gradient making the transition 
energy the same for all nodes. 

There are two important decoherence mechanisms which 
are given by the spontaneous decay of the involved Rydberg 
states and by the atomic motion. Exemplarily, for the 7 IS and 
61D states of rubidium and a distance of 20/im between nodes 
we calculate a transfer time of r =145 ns between two neigh- 
boring sites, radiative lifetimes including black-body radiation 
of > 100 /is and residual thermal motion that leads to a change 
of the interatomic distance of 1.4 fim per 100 /is at a tempera- 
ture of 1 /iK. Another source of decoherence is the interaction- 
induced motion ll22ll . We can model this motion quantita- 
tively |23] and calculate negligible changes of the interatomic 
distances of less than 0.2 /im per 100 /is. This means that both 
the chain and the elementary atomic system sustain coherence 
over timescales on the order of several ten /is and longer 

In conclusion, we have identified different time domains in 
the CTQW exciton decay in the presence of traps, domains 
which are directly related to the complex spectrum of the sys- 
tem's Hamiltonian. The CTQW average survival probability 
Hm {t) for an exciton to stay inside a linear system of N nodes 
with traps at each end can clearly be distinguished from its 
classical CTRW counterpart, PM{t). Finally, we proposed an 
experimental test for coherence on the basis of ultra-cold Ryd- 
berg atoms. 
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